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Abstract. The most interesting current open question in the theory of GRB afterglow 
is the propagation of jetted afterglows during the sideway expansion phase. Recent 
^^ , numerical simulations show hydrodynamic behavior that differs from the one suggested 

^^ ' by simple analytic models. Still, somewhat surprisingly, the calculated light curves show 

a 'jet break' at about the expected time. These results suggest that the expected rate 
of orphan optical afterglows should be smaller than previously estimated. 
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1 Introduction 

Our understanding of GRBs has been revolutionized by the BeppoSAX discovery 
of GRB afterglow. While GRBs last seconds or minutes the afterglow lasts days, 
weeks months or even years. This makes afterglow observations much richer. 
^^ , These observations provide us with multi-wavelength and multi-timescales data. 

t~^ ' At the same time the afterglow, which is a blast wave propagating into the 

^^ , surrounding matter is a much simpler phenomena than the GRB and it is possible 

to construct a simple theory that can be compared directly with the observations. 
-^ , In this short review we describe the theory of GRB afterglow. We begin with 

1-5 ' the simplest idealized model and continue with various levels of complications. 

The final level is full numerical simulations. We present preliminary results of 
O ' such simulations and compare them with analytic models. At present there is no 
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simple analytic explanation for the features seen in the numerical results. 



2 Spherical Hydrodynamics 



^^ ' The theory of relativistic blast waves has been worked out in a classical paper 

H i by Blandford & McKee (BM) already in 1976 0. The BM model is a self-similar 

spherical solution describing an adiabatic ultra relativistic blast wave in the limit 
r ^ 1. The basic solution is a blast wave propagating into a constant density 
medium. However, Blandford and McKee also describe in the same paper a 
generalization for varying ambient mass density, p = AR~'', R being the distance 
from the center. The latter case would be particularly relevant for fc = 2, as 
expected in the case of wind from a progenitor, prior to the GRB explosion. 

The BM solution describes a narrow shell of width ~ R/T^, in which the 
shocked material is concentrated, where F is the typical Lorentz factor. The 
conditions in this shell can be approximated if we assume that the shell is ho- 
mogeneous. Then the adiabatic energy conservation yields: 

E = JL.AR^-T^c^ , (1) 

O rv 
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where E is the energy of the blast wave and J7 is the soUd angle of the afterglow. 
For a full sphere f2 = 47r, but it can be smaller if the expansion is conical with 
an opening angle 9: fi = 2tt9^ (assuming a double sided jet). 

A natural length scale, / = [(3 — k)E/ilAc'^'\ , appears in equation |l[ 

For a spherical blast wave [2 does not change with time, and when the blast 
wave reaches R = I it collects ambient rest mass that equals its initial energy, 
the Lorentz factor F drops to 1 and the blast wave becomes Newtonian. The 
BM solution is self-similar and assumes F ^ 1. Obviously, it breaks down when 
R ~ I. We therefore expect that a Relativistic-Newtonian transition should take 
place around ^nr — l/c ~ 1.2yr(£^iso,52/'^i)^^^, where the scaling is for fc = 0, 
£^52 is the isotropic equivalent energy, Eiso — AttE/ fi, in units of lO^^ergs and 
ni is the external density in cm~^. After this transition the solution will turn 
into the Newtonian Sedov- Tailor solution. Clearly this produces an achromatic 
break in the light curve. 

The adiabatic approximation is valid for most of the duration of the fireball. 
However, during the first hour or so (or even for the first day, for fc = 2), the 
system could be radiative (provided that ee ~ 1). During a radiative phase the 
evolution can be approximated as: 

E = ^^AR^-'^FFoc^ , (2) 

where Fq is the initial Lorentz factor. Cohen, Piran & Sari derived an analytic 
self-similar solution describing this phase. Cohen & Piran g describe a solution 
for the case when energy is continuously added to the blast wave by the central 
engine, even during the afterglow phase. A self-similar solution arises if the 
additional energy deposition behaves like a power law. This would arise naturally 
in some models, e.g. in the pulsar like model M. 

3 Spherical Afterglow Models 

A good model for the observed emission from spherical blast waves can be ob- 
tained by adding synchrotron radiation to these hydrodynamic models. Sari, 
Piran & Narayan [^ used the simple adiabatic scaling (^ together with syn- 
chrotron radiation model and the relation between the observer time t, and R: 

t = R/CicF^ , (3) 

where Ci is a constant that may vary from 2 to 16 |q]. 

Assuming a powerlaw energy distribution of the shocked relativistic electrons: 
N{Eg) (X E~P, and that the electrons and the magnetic field energy densities 
are Ce and €b times the total energy density, Sari, Piran & Narayan B estimate 
the observed emission as a series of power law segments (PLSs), where 

F, oc t-'^iy-^ , (4) 

that are separated by break frequencies, across which the exponents of these 
power laws change: the cooling frequency, Vci the typical synchrotron frequency 
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Vm and the self absorption frequency Vsa- The analytic calculations were done 
for a homogeneous shell and for emission from a single representative point. At 
a specific frequency one will observe a break in the light curve when one of these 
break frequencies passes the observed frequency. An intriguing feature of this 
model is that for a given PLS, say for emission above the cooling frequency, 
there is a unique relation between a, P and p. The power law index p is expected 
to be a universal quantity as it depends on the, presumably common, acceleration 
processes and it is expected to be between 2 and 2.5 Ql- The consistency of those 
observed parameters could be a simple check of the theory. 

The simple solution, that is based on a homogeneous shell approximation, 
can modified by using the full BM solution and integrating over the entire vol- 
ume of shocked fluid [^. Such an integration can be done only numerically. It 
yields a smoother spectrum and light curve near the break frequencies, but the 
asymptotic slopes, away from the break frequencies and the transition times, 
remain the same as in the simpler theory. 

Chevalier & Lee |^ estimated the emission from a blast wave propagating into 
a wind profile n{R) ex R~^. They use equation (m) and calculate the synchrotron 
emission from a single representative point. This leads to different temporal 
scalings a of the PLSs, while the spectral indices /3 remain the same, since they 
are independent of the hydrodynamic solution. This results in different relations 
between a, (3 and p, providing in principle a way to distinguish between different 
neighborhoods of GRBs and between different progenitor models. 

Another modification to the " standard" model arises from a variation of the 
emission process. Sari & Esin flQ] considered the influence of Inverse Compton 
on the observed spectrum. They find that in some cases the additional cooling 
channel might have a significant effect on the observed spectrum and light curves. 



4 Jets 

The afterglow theory becomes much more complicated if the relativistic ejecta is 
not spherical. To model jetted afterglows we consider relativistic matter ejected 
into a cone of opening angle 9. Initially, as long as F ^ 9^^ |ll| the motion 
would be almost conical. There isn't enough time, in the blast wave's rest frame, 
for the matter to be affected by the non spherical geometry, and the blast wave 
will behave as if it was a part of a sphere. When F = C29^^ , namely atPl: 

1/3 / /) \ 8/3 



1 


/A 


( ^ ^ 


2(4-fc) 
(3-fc) 


Iday 


/^iso,52 


Ci 


\cj 


\C2J 


1 


ccl" 


V ni 



rapid sideway propagation begins. The last equality holds, of course for fc = 0. 

The sideways expansion continues with 9 ^ F^^. Plugging this relations 
in equation P) we find that R ~ const. This is obviously impossible. A more 

1 The exact values of the uncertain constants C2 and Ci are extremely important as 
they determine the jet opening angle (and hence the total energy of the GRB) from 
the observed breaks, interpreted as fjct, in the afterglow light curves. 



Tsvi Piran and Jonathan Granot 



detailed analysis il2,p3 14 1 reveals that according to the simple one dimensional 
analytic models F decreases exponentially with i? on a very short length scale n 

The sideways expansion causes a change in the hydrodynamic behavior and 
hence a break in the light curve. Additionally, when F ^ 9^^ relativistic beaming 
of light will become less effective. This would cause an extra spreading of the 
emission (that was previously focused into a narrow angle 9 and is now focused 
into a larger cone of opening angle F~^). If the sideways expansion is at the 
speed of light than both transitions would take place at the same time [^ . If 
the sideways expansion is at the sound speed then the beaming transition would 
take place first and only later the hydrodynamic transition would occur ||l(| . 
This would cause a slower and wider transition with two distinct breaks, the 
first and steeper break when the edge of the jet becomes visible and later a 
shallower break when sideways expansion becomes important. 

The analytic or semi-analytic calculations of synchrotron radiation from jet- 
ted afterglows [[l2|jl5| , pij| , [l7| , |l4[ have led to different estimates of the jet break 
time ijet and of the duration of the transition. Rhoads iQ calculated the light 
curves assuming emission from one representative point, and obtained a smooth 
'jet break', extending ~ 3 — 4 decades in time, after which -Fl/>i/„ ex t~^ . Sari 
Piran & Halpern |g5[ assume that the sideway expansion is at the speed of light, 
and not at the speed of sound (c/a/S) as others assume, and find a smaller 
value for tjct- Panaitescu and Meszaros p6| included the effects of geometrical 
curvature and finite width of the emitting shell, along with electron cooling, 
and obtained a relatively sharp break, extending ^1 — 2 decades in time, in 
the optical light curve. Moderski, Sikora and Bulik [^ used a slightly different 
dynamical model, and a different formalism for the evolution of the electron 
distribution, and obtained that the change in the temporal index a {F^ oc i~") 
across the break is smaller than in analytic estimates (a = 2 after the break for 
u > Vm, p = 2.4), while the break extends over two decades in time. Kumar and 
Panaitescu El find that for a homogeneous (or stellar wind) environment there 
is a steepening of Aa ~ 0.7 (0.4) when the edge of the jet becomes visible, while 
the steepening due to sideways expansion extends over 2 (4) decades in time. 
They conclude that a jet running into a stellar wind will not leave a prominent 
detectable signature in the light curve. 

The different analytic or semi-analytic models have different predictions for 
the sharpness of the 'jet break', the change in the temporal decay index a across 
the break and its asymptotic value after the break, or even the very existence 
a 'jet break' [|l^. All these models rely on some common basic assumptions, 
which have a significant effect on the dynamics of the jet: (i) the shocked matter 
is homogeneous (ii) the shock front is spherical (within a finite opening angle) 
even at t > tj^t (iii) the velocity vector is almost radial even after the jet break. 

However, recent 2D hydrodynamic simulations [|l9| show that these assump- 
tions are not a good approximation of a realistic jet. Figure shows the jet at 

^ Note that the exponential behavior is obtained after converting equation ry to a dif- 
ferential equation and integrating over it. Different approximations used in deriving 
the differential equation lead to slightly different exponential behavior, see ]l3|. 
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Fig. 1. A relativistic jet at the last time step of the simulation |1!J|. (left) A 3D view of 
the jet. The outer surface represents the shock front while the two inner faces show the 
proper number density {lower face) and proper emissivity {upper face) in a logarithmic 
color scale, (right) A 2D 'slice' along the jet axis, showing the velocity field on top of 
a linear color-map of the lab frame density. 



the last time step of the simulation. The matter at the sides of the jet is prop- 
agating sideways (rather than in the radial direction) and is slower and much 
less luminous compared to the front of the jet. The shock front is egg-shaped, 
and quite far from being spherical. Figure shows the radius R, Lorentz factor 
r, and opening angle of the jet, as a function of the lab frame time. The rate 
of increase of 9 with R « ciiab, is much lower than the exponential behavior 
predicted by simple models |E2| . The value of 9 averaged over the emissivity is 
practically constant, and most of the radiation is emitted within the initial open- 
ing angle of the jet. The radius R weighed over the emissivity is very close to 
the maximal value of R within the jet, indicating that most of the emission orig- 
inates at the front of the jein, where the radius is largest, while R averaged over 
the density is significantly lower, indicating that a large fraction of the shocked 
matter resides at the sides of the jet, where the radius is smaller. The Lorentz 
factor r averaged over the emissivity is close to its maximal value, (again since 
most of the emission occurs near the jet axis where F is the largest) while F 
averaged over the density is significantly lower, since the matter at the sides of 
the jet has a much lower F than at the front of the jet. The large differences 
between the assumptions of simple dynamical models of a jet and the results of 
2D simulations, suggest that great care should be taken when using these models 
for predicting the light curves of jetted afterglows. Since the light curves depend 
strongly on the hydrodynamics of the jet, it is very important to use a realistic 
hydrodynamic model when calculating the light curves. 

Granot et al. p9[ used 2D numerical simulations of a jet running into a con- 
stant density medium to calculate the resulting light curves, taking into account 
the emission from the volume of the shocked fluid with the appropriate time de- 
lay in the arrival of photons to different observers. They obtained an achromatic 
jet break for v > i^m{t jet) (which typically includes the optical and near IR), 

^ This implies that the expected rate of orphan optical afterglows should be smaller 
than estimated assuming significant sideways expansion! 
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Fig. 2. The radius R {left frame), Lorentz factor F — 1 (middle fram,e) and opening 
angle 6 of the jet {right frame), as a function of the lab frame time in days |19| ]. 

while at lower frequencies (which typically include the radio) there is a more 
moderate and gradual increase in the temporal index a at tjet, and a much more 
prominent steepening in the light curve at a latter time when v^ sweeps past 
the observed frequency. The jet break appears sharper and occurs at a slightly 
earlier time for an observer along the jet axis, compared to an observer off the 
jet axis (but within the initial opening angle of the jet). The value of a after the 
jet break, for v > v„i, is found to be slightly larger than p {a = 2.85 for p = 2.5). 
Somewhat surprisingly we find that in spite of the different hydrodynamic 
behavior the numerical simulations show a jet break at roughly the same time 
as the analytic estimates. This encourages us to trust the current estimates of 
the jet opening angles. However, we should search for an intuitive explanation 
for the nature of the hydrodynamic behavior and for a simple analytic model 
that would predict it. 
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